Vertical probe-induced asymmetric dust oscillation in complex plasma 
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Spherical, /im-sized particles within a Coulomb crystal levitated in the sheath above the powered 
lower electrode in a GEC reference cell are perturbed using a Zyvex S100 Nanomanipulator. Using 
the S100, a vertical probe is positioned within the cell at various locations with respect to the 
crystal formed within the sheath. As the probe is lowered toward the horizontal plane of the dust 
layer, a circular cavity opens in the center of the crystal and expands. To explore the minimally 
perturbative state, the probe is lifted to the position that closes this cavity, the probe potential is 
oscillated, and the motion of the particle directly beneath the probe is analyzed. Using a simple 
electric field model for the plasma sheath, the change predicted in the levitation height is compared 
with experiment. 

PACS numbers: 52.27.Lw 



I. INTRODUCTION 

When macroscopic dust is introduced into a plasma, 
the resulting complex plasma system offers the oppor- 
tunity to simulate, at its most fundamental level, in- 
teresting physics pertaining to a host of research areas, 
including protoplanetary [1 or protostellar [2 develop- 
ment, contamination in plasma enhanced semiconductor 
atomic layer deposition and etching systems [3] , and wall 
erosion within a fusion device [4]. At low temperature 
and power, it also provides a basis for the study of en- 
tirely new plasma physics. 

A recurring problem in complex plasma physics is 
obtaining accurate measurements in the sheath region, 
which is nonuniform and spatially limited compared to 
the bulk. Dust particles in a capacitively coupled rf 
plasma levitating beneath the bulk plasma can act as in 
situ probes providing insight into the fundamental sheath 
parameters [5]. One difficulty with this method is accu- 
rately determining the charge on the dust particles, which 
is itself linked to the plasma parameters. Given known 
plasma conditions, Orbital Motion Limited (OML) the- 
ory [6] can be used to predict the amount of charge col- 
lected by each particle. An improved charge estimation 
can be obtained using a shifted Maxwellian distribution 
to incorporate the effect of ions streaming toward the 
lower electrode [7 , although this requires an estimation 
of the ion velocity which is also not easily measured. 

A plate with a milled depression 1-mm deep placed on 
the lower electrode can provide horizontal confinement, 
allowing dust particles to be stably levitated in a GEC 
rf cell through a balance of the vertical forces. In this 
case, the dust can achieve a minimum energy state by 
organizing into a two-dimensional circular lattice with 
hexagonal symmetry. 
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FIG. 1. Dust particle positions in a planar dust crystal with 
an open central cavity created by a negatively biased vertical 
probe located at the center of the cavity region. This was 
achieved using a system power of 10 W, a neutral gas pres- 
sure of 100 mTorr, a probe potential of -30 V, and a 1-inch 
diameter circular cutout for horizontal confinement. 



Many methods have been used to manipulate the dust 
in such a complex plasma crystal. These include modi- 
fying the dc bias on the lower electrode [8 , applying ra- 
diation pressure [9], initiating a wave through the crys- 
tal using a powered horizontal wire [10], instituting a 
temperature gradient to establish a thermophoretic force 
[TT] . altering the particle charge or sheath electric field 
by changing basic plasma parameters such as the power 
and/or neutral gas pressure [12] . and beaming a pulsed 
stream of ions [13] . 

The addition of a probe to the system provides the abil- 
ity to apply a controlled potential difference at various 
positions relative to the charged dust layers. Previous 
experiments have used such probes to induce voids in a 
DC plasma [14] and in microgravity rf plasmas [T5] . 
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FIG. 2. (Color online) The Zyvex S100 Nanomanipulator 
head is located inside the plasma cell, mounted upside down 
above the upper ring electrode. One of the four manipulators 
shown is used in this experiment to position the probe as 
described in the text. 



Parameter 


Settings 


Probe Bias 
Probe Peak to Peak 
System Power 

Pressure 
Probe Height 
Frequency 
DC Bias 


10, 20, 30 V 
45, 55, 65 V 
0.75, 1.00, 1.50 W 
70, 80, 90 mTorr 
7, 9, 11 mm 
1, 2.3, 5 Hz 
-10, -5, -1 V 



TABLE I. Experimental parameters used in this experiment. 
The bolded value for each parameter defines the base config- 
uration, as discussed in the text. 



In this experiment, an adjustable probe is positioned 
perpendicularly to a horizontal dust lattice. As it ap- 
proaches the crystal plane, the particles move radially 
outward, creating an open circular space within the crys- 
tal as shown in Fig. [T] This will henceforth be referred 
to as the cavity. Under operating conditions that just 
close the cavity, a numerical model is used to calculate 
the charge on the dust particle. The probe potential is 
oscillated, and the resultant dust motion is analyzed to 
determine the neutral drag coefficient and the resonant 
frequency. The plasma parameters within the model are 
then varied to find the best fit to observed particle os- 
cillations. Assuming probe interaction with the dust is 
a superposition of near (i.e., direct interaction with the 
probe potential) and far (i.e., changes to the plasma) field 
effects, this experiment isolates the latter. 

The arrangement of this paper is as follows: First, the 
experimental apparatus is described. Second, the exper- 
iment itself is introduced and forced particle oscillation 
over a range of parameters is explored. Third, a simple 
electric field model is applied to predict the equilibrium 
levitation heights of the crystals. Different parameters 
of this model are then perturbed to incorporate the os- 
cillation, and the results are compared to experiment. 




FIG. 3. Side view of the vertical oscillation experiment. For 
the 'base configuration' setup, the tip of the probe is located 
7.3 mm above the dust equilibrium levitation height. Only the 
particle between the (dashed) bars is tracked, with the area in 
the box, which excludes the electrodes, used for calibration. 



Finally, the conclusions are reported. 



II. 



APPARATUS 



A modified GEC (Gaseous Electronics Conference) 
rf reference cell [T6l [17] located at the Center for As- 
trophysics, Space Physics, and Engineering Research 
(CASPER) was employed for this study. The cell con- 
tains two electrodes, 8 cm in diameter, separated by 
1.9 cm. The lower electrode is powered at 13.56 MHz 
while the upper ring-shaped electrode and chamber are 
grounded. In this experiment 8.9 micron diameter 
melamine formaldehyde (MF) dust particles were levi- 
tated in an argon plasma. A plate containing a circular 
cutout 1 mm deep and 25.4 mm (1 inch) in diameter 
was placed on the lower electrode to provide horizontal 
confinement of the dust particles. Images of the dust 
particles were captured at 125 and 250 frames per sec- 
ond using a side-mounted CCD camera. Table |l] lists the 
system parameters used for this experiment. The actual 
power delivered to the plasma, calculated by measuring 
the phase difference between the current and the volt- 
age, was 0.23, 0.36, and 0.56 W (U r/ =9.6, 11.4, and 14.1 
V). The DC bias of the lower electrode was held fixed 
using an external power supply. Langmuir probe mea- 
surements of the bulk plasma were also collected using a 
SmartProbe (Scientific Systems Ltd) [T7] . 

The Zyvex S100 [18] head (originally designed for sam- 
ple manipulation in microscopy applications) was at- 
tached within the plasma chamber to act as a pertur- 
bation tool (Fig. [2]). A hollow cylindrical probe 48 mm 
in length, having an outer diameter of 450 /im, was con- 
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FIG. 4. (Color online) Normalized optical emission intensity 
profile (solid blue line) of the plasma (Fig. |3| averaged over 
the horizontal coordinate as a function of the distance above 
the lower electrode (where the upper electrode is located at 
z=19 mm). The derivative of the intensity (green dots) is 
superimposed with vertical lines indicating the position of the 
dust (red dashed), the 1/e point (green dash-dot), the local 
maximum of the derivative emission intensity (dotted black), 
and the maximum of the intensity (solid blue). 
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FIG. 5. (Color online) Distance from the lower electrode ver- 
sus time for: the particle height (red circles), the sheath edge 
(green triangles) , the maximum of the derivative of the emis- 
sion intensity (black line), and the maximum of the emission 
intensity (blue crosses), for system powers of a) 0.75 W, b) 
1.00 W, and c) 1.50 W. 



nected to the S100 head. A tip, with an adjustable length 
of up to 14 mm and a diameter tapering from 250 to 50 
microns over the last 100 microns, protruded from the 
probe (Fig. [3]). The S100 is capable of remote controlled 
movement of the probe by up to 10 mm in all three di- 
mensions. The probe potential can be controlled with 
respect to the ground using an external power supply. 



III. EXPERIMENT 

Motivation for the current experiment stemmed from 
the observation that changing the probe bias and height 
both produced and altered the size of a circular cavity 
in the dust layer. In the horizontal wire experiment ref- 
erenced above [10] , the primary repulsion force upon the 
dust was shown to be produced by the direct electric 
field. This proves not to be the case here. In all cases 
the cavity created is much larger than would be expected 
assuming a simple interaction between the probe and the 
dust through a shielded Coulomb potential. In order to 
examine this behavior, the probe tip was positioned at 
the transition height, defined as the point of closest ap- 
proach before a cavity is opened in the crystal. A side 
view showing the probe at the transition height for a 
crystal at its equilibrium levitation height for the base 
configuration as defined in Table [I] is shown in Fig. [3J 

Datasets consisting of 125 images taken at 125 fps were 
obtained for each of the parameters shown in Table [l| 
The set of parameters in bold (20 V probe bias, 55 V 
probe potential peak to peak, 1.0 W system power, 80 



mTorr gas pressure, 7.3 mm probe height, 2.3 Hz oscilla- 
tion frequency, and -5 V fixed DC bias) will henceforth 
be referred to as the base configuration. As each param- 
eter is modified in the experiment, all others remain the 
same as those established for the base configuration. 

Images were analyzed employing both particle tracking 
[T9] and profile analysis [20] . Image calibration was com- 
pleted by subtracting the minimum intensity of the image 
set in the region excluding the electrodes (indicated by 
the box in Fig. |3| from all pixels before determining the 
position of the maximum intensity and sheath edge. The 
probe was removed from the image before profile analysis 
by copying the adjacent region. 

The quantities best characterizing the oscillation are 
indicated in a representative image profile of the base 
configuration in Fig. [4] These include, in ascending ver- 
tical position, the particle levitation height, the sheath 
edge (defined as the position where the optical emission 
intensity decreases by a factor of e from its maximum 
[21]), the location of the maximum derivative of the emis- 
sion intensity in the lower sheath, and the point of emis- 
sion maximum. Sinusoidal oscillation of the probe po- 
tential was found to yield non-sinusoidal change in these 
four values, as shown in Fig.[5j over two oscillation cycles 
for three different system powers. 

The fact that the maximum emission intensity moves 
lower with an increase in power is most likely due to 
cell geometry; the plasma exhibits an asymmetric vertical 
emission profile (averaged over the horizontal direction) 
as shown in Fig. [4] due in part to the fact that the upper 
electrode is a hollow ring and grounded whereas the lower 
electrode is a powered plate. Oscillation of the probe 
potential causes the change in amplitude of the emission 
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FIG. 6. (Color online) Particle oscillation amplitude versus 
probe oscillation frequency. The fit shown is for a damped, 
forced harmonic oscillator. 



intensity maximum to become significant at lower powers 
([5^l) , but reduced at higher powers (J^). 

It was found that the position of the local maximum of 
the derivative of the emission profile of the lower sheath 
remains constant with change in power (Fig. [5]). This is 
true not only for variations in power, but over all param- 
eters tested. As shown in Fig. [4j the local maximum is 
located at z=6.6 mm, with a local minimum at z=l6.8 
mm in the upper sheath. The peak which can be seen at 
z=4.7 mm is due to the line of dust. 

The location of the 1/e point, indicating the sheath 
edge, decreases with an increase in power while the par- 
ticle levitation height increases (Fig. [5J, both due to 
an increase in ionization of the plasma. The defini- 
tion of sheath edge position used in this experiment was 
first proposed by Beckers et al. [21]. An alternative ex- 
perimental method to determine the sheath edge based 
on theory is measurement of the equilibrium height of 
nanoparticles [22 . Applying the definition of the sheath 
edge as the 1/e point to Fig. 4 in [22] results in 0.7% 
agreement in vertical position, justifying the choice in 
the present work. 

The response of the dust to the probe oscillation can be 
approximated as in [8] by the steady-state response of a 
damped, forced harmonic oscillator, with the amplitude 
given by 

F 

A ^ = [(^-^)2 +4/ 32 CJ 2]l/2 W 

where f3 is the damping coefficient, ujq is the resonant 
frequency, uj is the driving frequency, and F is the mag- 
nitude of the driving force divided by the mass of the 
dust particle. A plot of the dust particle amplitude ver- 
sus frequency (Fig. [6| shows good agreement with theory. 
At a pressure of 80 mTorr, the resulting fit parameters 
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FIG. 7. (Color online) Expanded view of the base configu- 
ration (Fig. [5)3). The solid (red) line is a polynomial fit to 
the particle position (red circles) over time while the dashed 
line represents the fit to the sheath position (green squares). 
The vertical lines shown indicate the peaks of these fit lines il- 
lustrating the delay between maximum particle response and 
maximum sheath response. For this case, the delay is 15.0 ms; 
delays of similar magnitude were found for other experimental 
parameters. All values are shown as percent differences from 
their respective minima. 

are /3=9.66 s -1 , cj =65.0 rad/s (10.3 Hz), and F=1.08 
N/kg. These are comparable to the results reported in 
[23] . where f3=8 s _1 , and cjo/27t=13.2 Hz, measured at 66 
mTorr. The process that creates the non-sinusoidal par- 
ticle response in time must not be sensitive to changes 
to frequencies in this range, and thus the only impact to 
Eqn. [I] occurs in F. 

The phase delay which can be seen between the max- 
ima of the sheath edge and the particle positions as shown 
in Fig. [7] is well known for damped harmonic oscillation. 
For all cases, it is found to be in the correct direction with 
the particle response lagging the driving force. While the 
probe potential oscillates sinusoidally, the location of the 
sheath edge does not, as seen in Fig. [5] This lends a hy- 
pothesis that the force on the particle is not sinusoidal. 
Though a periodic, non-sinusoidal driving force will yield 
the same resonant frequency as that derived using Eqn. [I] 
[T2] , calculation of the phase shifts under different driving 
forces requires a numerical approach. The phase differ- 
ence (S) for a sinusoidal force, 

applied to this experiment does not lead to a damping 
constant (12.7 s _1 ) or resonant frequency (7.7 Hz) con- 
sistent with that found from the frequency sweep (Fig. [6]). 
Deviation of the experimental data from the fit line at 
half the resonant frequency shown in Fig. [6] was also seen 
in two other experiments. Ivlev et al. [24] drove verti- 
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FIG. 8. (Color online) Intensity space-time contour plot for 
the base configuration. The colorbar gives the intensity rela- 
tive to the overall minimum, and normalized to the resulting 
maximum ((/ — I m in) / Imax)- The sinusoidal probe potential 
is superimposed. 



cal dust oscillations, using a powered wire placed below 
the particle layer and oriented horizontally (i.e., parallel) 
to the lower electrode. While the amplitude fit deviated 
from a damped harmonic oscillator above driving poten- 
tials of 50 mV peak to peak, notable superharmonic re- 
sponse appeared at driving potentials (4 V peak to peak) 
substantially higher than those used in this experiment. 
Homann et al. [25] examined particle response to lower 
electrode DC bias oscillations of square and sine waves. A 
laser was also pulsed on the particle without bias changes 
which created a square wave force without sheath mod- 
ification. Superharmonic response was found only when 
the force delivered to the dust was not sinusoidal. There- 
fore the phase shift discrepancy and superharmonic peak 
are two results that confirm that the force on the par- 
ticle is not sinusoidal and indicate that changes in the 
plasma are driving the particle's motion, rather than the 
sinusoidal potential on the probe. 

Changes in the plasma can be characterized through 
the plasma intensity. Fig. [8] shows the evolution of the 
plasma intensity over time. The vertical position where 
the change in intensity is most pronounced is in the 
plasma bulk. A decrease in the plasma glow results from 
a reduction in electron density, which in turn results in a 
reduction of the total argon electron transitions. When 
the probe reaches its maximum positive potential, the 
bulk plasma intensity decreases by up to 9%, as shown 
in Fig. [9] The glow responds faster than the particle, 
providing an indicator for the force on the particle. This 
presents a link to the particle delay (Fig. [7]) without hav- 
ing to independently match the probe potential to the 
oscillation. Fixed positive probe potentials corresponded 
to the dust crystal height being raised as a whole, and the 
height of the maximum optical emission intensity from 



FIG. 9. (Color online) Maximum intensity percentage change 
over time for the base configuration (solid blue line), super- 
imposed with tracked particle data (dashed green line) and 
probe potential (dash-dot black line), which is scaled to the 
maximum particle position. 
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FIG. 10. (Color online) Oscillation amplitudes for the sheath 
edge (green squares) and the particle (red circles) with 
quadratic or linear fits for variation in individual parameters 
given in Table]!] a) Probe bias (V), b) probe peak to peak (V), 
c) system power (W), d) pressure (mTorr), e) probe height 
(mm), f) frequency (Hz), g) DC bias (V). 



the plasma shifted upward as well. 

The dust oscillation amplitudes for the parameters 
tested is shown in Fig. [To] The amplitude of the sheath 
oscillation shown is larger than that observed for the par- 
ticle oscillation. This implies that the particle is not sim- 
ply entrained in the plasma, but instead that its reduced 
amplitude is due to its mass (resisting sudden change 
through inertia); cross-sectional area (motion impeded by 
colliding neutral gas atoms and ions); charge (Coulomb 
collisions with the streaming ions); and the structure of 
the sheath electric field. While any increase in probe 
bias generates a nonlinear increase in oscillation ampli- 
tude (lOi), an increase in the peak to peak oscillation 

(flOlb). On the 



voltage results in a nearly linear increase 
other hand, increasing the rf power causes a decrease in 
the oscillation amplitude (10 z) due to ionization increase 



and electron temperature decrease (as found by Lang- 
muir probe measurements), both of which increase the 
overall shielding by decreasing the Debye length. There- 
fore, a reduced perturbation is applied to the plasma, 
resulting in a smaller particle amplitude. At the same 
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time, changing the pressure and probe height do not ap- 
preciably affect the oscillation amplitude (101 and [lOfc) . 
As the probe oscillation frequency approaches 5 Hz, the 
particle no longer has time for relaxation between cy- 
cles and the maximum particle amplitude increases due 
to superposition ([Tof ). Decreasing the DC bias on the 
lower electrode (|1(%) increases the potential difference 
between the plasma and the electrode, raises the parti- 
cle within the sheath region, and results in a nonlinear 
increase in the particle amplitude. Taken together, these 
results indicate once again that the probe locally modi- 
fies the plasma which in turn drives the oscillation of the 
particle, as supported by the phase delay shown in Fig. [7] 



IV. 



MODEL 



A numerical model was developed to examine the data 
described above. The general equation of motion for a 
dust particle is given by 



mz = F E - F q - Fi + Fr 



(3) 



where F n are the forces produced by the electric field, 
gravity, streaming ions, and neutral damping (whose sign 
changes so that it always opposes the particle motion), 
respectively. The electric potential within the sheath is 
assumed to be parabolic [26 , where the fixed DC po- 
tential on the lower electrode is Vq at z=0 and equal to 
the plasma potential (V p ) at the sheath edge (z=d). The 
electric force acting on a dust particle is then: 



Fe{z) = 



2Q(V P - Vo) 
d 



(s- 1 ) 



(4) 



The charge (Q) on each particle is initially assumed to 
be 14,750e, where e is the elementary electron charge, as 
determined by previous experiment in this GEC cell [23] . 
This model was validated by first using it to predict the 
dust equilibrium levitation height, and then to estimate 
the charge on a dust particle for each configuration. Af- 
ter the charge values were found, changes in the plasma 
based on changes of the experimental parameters were 
explored. 

It is well known that ions stream from the plasma bulk 
and are accelerated into the sheath [27]. For this exper- 
iment, it is assumed these ions exit the sheath edge at 
the Bohm velocity into a collisionless sheath. However, 
the ion-neutral mean free path (625 microns at 80 mTorr 
using the total cross-section from [28 ) is smaller than 
the sheath thickness, which means the sheath exhibits 
some collisionality. Ignoring this effect overestimates the 
ion drag on a dust particle, since some momentum is 
transferred to the neutrals; however, this is still quite 
small relative to the other forces considered [21 . Ion 
drag was calculated using the model of Khrapak et al. [7] 
and added to the downward gravitational force. The ion 
drag calculation was modified to determine the ion and 
electron densities at the dust position given the electric 
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FIG. 11. (Color online) Forces involved in the levitation of 
the dust within the plasma sheath. The dashed (red) hori- 
zontal line represents the downward force of gravity including 
ion drag (hence the small deviation from -1). Using the ini- 
tial estimate for the particle change in Eqn.[4]gives the linear 
electric field force in the sheath (blue solid diagonal line). 
The intersection of these two lines gives the predicted equi- 
librium height (black solid vertical line). The experimentally 
measured levitation height (black dashed line) can be used to 
determine the actual particle charge by adjusting the electric 
force (blue dash-dot line) so that it intersects the combined 
downward force at the levitation height. 
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FIG. 12. (Color online) Percent difference of the predicted 
height from the experimental equilibrium position, based on 
the assumed charge for the base configuration. 



potential at the dust height, using the densities in the 
plasma bulk. These modifications included energy con- 
servation and the continuity equation for the ions [29] 
(which quantifies the decrease in ion density as the ions 
accelerate toward the lower electrode), and a Boltzmann 
distribution of electrons. 

The initial force balance for the base configuration (as 
given in Table |l| is shown in Fig. 11 , where the forces 
are normalized by the gravitational force. Application 



of the model under the parameter variation used is seen 
in Fig. [l2j where the predicted height (using Q=14,750e 
for each case) differs from the actual levitation height 
by less than 3%. The charge on the particle was then 
adjusted so that the electric force balances the grav- 
itational and ion forces at the experimental levitation 
height. This is 16,590e for the base configuration, which 
is 12% larger than the initially assumed charge (generat- 
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Parameter 


Low Mid High 


Parameter 


Low Mid High 


Probe Bias 
Peak to Peak 
Power 
Pressure 


14.9 16.6 18.2 
16.1 16.6 16.1 
10.7 16.6 39.3 
18.0 16.6 15.8 


Height 
Frequency 
DC Bias 


16.6 15.3 16.9 

16.4 16.6 16.2 

14.5 16.6 18.2 



TABLE II. Dust particle charge found from the balance of the 
forces discussed in the text, in units of 10 3 electron charges. 
All parameters were changed independently relative to the 
base configuration. The large particle charge for the highest 
power occurs due to the particle approaching the sheath edge. 



ing a 2% change in height). Note that the experimental 
charge must be found numerically because changing the 
grain charge changes the ion drag force, which is given 
by a transcendental function, which does not admit an 
analytical solution. 

Table [TT] gives the resulting particle charges found for 
all parameters. As shown, the charge increases with 
power, DC bias, and probe bias, but decreases with pres- 
sure, consistent with the observed rise and fall in particle 
height, respectively. Calculating the percentage differ- 
ence in charge as the remaining parameters are changed 
finds them to be consistent within 8%. 

The nonlinear increase in charge found with increasing 
power may be due to overestimation of the ion speed as 
mentioned previously. In the P=1.5 W case, the particle 
levitates near the sheath edge. A nonlinearity in the 
electric field [30] occurs here which is not incorporated 
into the model. This results in an unrealistic increase 
in the ion drag. Note that the resulting charge increase 
with power is due to both the decrease in position of the 
sheath edge and the increase in the particle height as 
shown in Fig. [5] 

Comparing the experimental data with the model al- 
lows examination of the mechanism behind the asym- 
metric oscillation observed. There are several ways to 
perturb the electric field model in order to produce a 
vertical oscillation. First, the 1/e point, when experi- 
mentally tracked over time, can be used to adjust the 
position of the sheath edge, which changes the electric 
field calculated by Eqn. [4j Second, the potential at the 
sheath edge (i.e., the plasma potential) can be altered 
sinusoidally to follow the potential on the probe while 
the sheath edge remains fixed. Third, the charge of the 
particle can be changed by oscillating the plasma power, 
which also changes the ionization rate and plasma den- 
sity. Since the change in plasma power and the manner in 
which it affects the grain charge requires many assump- 
tions, the charge is instead estimated as a function of 
the grain's vertical position and the plasma emission in- 
tensity measured during probe potential oscillation. For 
each of these methods, the force balance can be used to 
predict the position of the particle over time. A change 
in the ion flow due to the probe potential was also consid- 
ered. However, when the probe is more positive than the 
plasma potential, downward ion flow should increase in 
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FIG. 13. (Color online) Predictions of multiple perturbations 
to the simple electric field model for the base configuration, a) 
The (black) dots give the the experimental particle position. 
The (red) squares represent the result of adjusting the sheath 
edge. The (blue) diamonds represent changing the plasma 
potential based on the probe oscillation. The (green) circles 
represent changing the grain charge (with the ion drag cal- 
culated from the equilibrium charge). The RMS deviation of 
the perturbations to the experimental data are shown in b)- 
h) where the parameter varied is b) probe bias (V), c) probe 
peak to peak (V) , d) system power ( W) , e) pressure (mTorr) , 
f) probe height (mm), g) frequency (Hz), and h) DC bias (V). 



speed, reducing the particle levitation height. Therefore, 
it can not independently generate this oscillation. 

The results of the perturbations discussed above are 
shown in Fig. [l3j using the experimentally determined 
particle charge. As shown, calculations based on apply- 
ing the experimental shift of the sheath edge yield the 
best prediction of the particle's motion (green circles). 
Taking the probe bias directly into account by applying 
a sinusoidal plasma potential results in sinusoidal parti- 
cle motion (red squares). For the base configuration, the 
probe bias and amplitude used in the model had to be 
reduced from experimental values by about 50% in order 
to match the measured equilibrium height and oscillation 
amplitude of the dust. This reduction is due to plasma 
shielding: a Debye length of 695 microns was calculated 
following the formula in [31] (which was used for the ion 
drag calculation). However, an increased shielding is ex- 
pected since the probe is several millimeters away from 
the dust. Because the electron temperature change in 
the vertical direction was not estimated (such as in [22 J, 
shielding was not included. For the third perturbation, 
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the motion is created by allowing the charge to change 
as a function of the maximum emission intensity, /(£), 
resulting in motion indicated by the blue diamonds. In 
this case, the charge ranges from its equilibrium value 
(Qo) to the value found from the force balance when the 
particle is at its maximum position (Q ma x)- This result 
is adjusted following an interpolation inversely propor- 
tional to the reduction in I(t) shown in Fig. [9j whose 
functional form is 



Q(J(t)) = Q + (Q 



QoV 



i(t) 



(5) 



where I min (I ma x) is the minimum (maximum) value over 
time of the maximum plasma intensity. The ion drag cal- 
culation in this case is fixed, and uses only the equilib- 
rium charge value. 

Although a reduction in intensity generally signals a 
reduction in power delivered to the plasma and thus a 
reduced grain charge, the fact that the particle position 
increases implies the reverse. There are several possi- 
ble explanations for this. Increasing the electron flow 
to the positive probe would increase the electron flux 
impacting the dust. The particle could also rise due 
to the increase in plasma shielding length due to a re- 
duced electron density; this would increase direct attrac- 
tion by the probe and repulsion from the lower electrode. 
The probe's electric field could interact with that of the 
sheath as well; however, since the tip is approximately 
10 shielding lengths away from the dust, these effects are 
limited. Therefore, the most likely process is that the 
plasma bulk contracts due to a reduction in electrons as 
the probe collects them, raising the sheath edge, while 
the charge remains relatively constant. 

The plots shown in Fig. 13 (b-h) illustrate the root 
mean square (RMS) deviations for the parameters tested 
as compared to experimental results. The RMS deviation 
for modifying the location of the sheath edge is smallest 
in 13 of the 15 cases, implying this is the best possible 
explanation. However, the other effects may contribute 
to the oscillation as well. Besides these perturbations, 
the speed of the flowing ions could be directly affected 
by the probe, but the ion drag would need to increase a 
few orders of magnitude in order to independently gen- 
erate the observed particle amplitude, and the greatest 
motion would then occur during the maximum negative 
potentials (a phase shift of tt from that seen in the ex- 
periment). Therefore this oscillation can be considered a 
far field effect, unlike the near field effect studied in the 
powered horizontal wire experiment [TO] . 



V. CONCLUSIONS 

An adjustable, powered vertical probe was used to pro- 
vide direct measurement of the neutral drag coefficient, 
the resonant frequency, and the charge on levitating par- 
ticles in an rf powered dusty plasma. Measured values 
are comparable to those found in previously published 



experiments [8] [23] . The probe potential was oscillated 
to create vertical oscillations of the dust particles. The 
amplitude of the resulting dust oscillation was shown to 
be most affected by increasing the probe bias; this re- 
sults in the probe potential spending more time above the 
plasma potential. Other parameters substantially alter- 
ing the dust oscillation amplitude were determined to be 
the probe's peak to peak oscillation potential, the lower 
electrode DC bias, and the plasma power (Fig. 10). 



A phase delay between the driving force and the maxi- 
mum particle height was observed (Fig. [7]). Such a phase 
delay is predicted for a forced, damped oscillator, with 
the calculation dependent on the forcing function. The 
measured particle oscillation amplitude was shown to fol- 
low that of a forced, damped harmonic oscillator, as seen 
by a fit to the amplitudes over a range of frequencies 
(Fig.|6|. Interestingly, the phase delay did not. Although 
a sine wave was applied to the probe, asymmetric oscil- 
lation of the particle and sheath edge showed the probe 
primarily interacts not directly with the particle but in- 
stead with the plasma. 

The plasma discharge intensity was shown to decrease 
by a maximum of 9% during probe oscillation (Fig. [9|, 
signaling a reduction in the electron density in the 
plasma. The greatest decrease occurs in the bulk (Fig.(8|. 
One unexpected result was that the position of the local 
maximum of the derivative of the emission profile in the 
lower sheath remained constant (Fig. [5| under variation 
in all system parameters. 

The electric field in the sheath was studied under ap- 
plied perturbation by analyzing the plasma glow. Un- 
der the low power regime examined, the 1/e point of 
the emission profile was shown to be a useful measure- 
ment of the sheath edge, as corroborated by [21] or con- 
firmed when applied to another experiment [22]. Cal- 
culating a force balance (Fig. 11) between the sheath 
electric field, gravity, and ion drag led to a predicted 
levitation height in agreement to that found in experi- 
ment (Fig. 12). Thus, the greatest change in dust parti- 
cle charge is caused by changes in plasma power, neutral 
gas pressure, and DC bias on the lower electrode. 

Perturbations to the electric field model examined 
three possible ways of generating this oscillation: making 
the sheath height dynamic, changing the plasma poten- 
tial, and adjusting the charge on the particle (Fig. [l3^). 
RMS calculations showed that changing the location of 
the sheath edge in time yielded the best fit to the col- 
lected particle data (Fig, plfo-h). 

Several areas of this experiment could be expanded. 
With appropriate calibration the intensity reduction ob- 
served in the bulk could be used to calculate the change 
in electron density. This could provide an extra parame- 
ter in plasma electromagnetic wave transmission studies, 
such as Faraday rotation. The fixed point of the position 
of the local maximum of the derivative of the emission 
profile may be purely dependent on the geometry of the 
cell and might be useful as a point of reference for image 
analysis. Theoretical work should be conducted to un- 
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cover what conditions can generate this result. Ion drag 
in the numerical model may be improved by changing 
the formulas used to calculate the decrease in ion and 
electron densities from the sheath edge [17], or testing 
the importance of the ion-neutral collisions and the non- 
linearity near the sheath edge. Though the ion drag is 
small for most of the plasma parameters in this experi- 
ment, it does increase significantly with power and would 



increase if the ion speed was found to be greater than the 
minimum Bohm velocity. 
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